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\q ! Abstract 
O 
t> 

' We find the leading order result for the anapole form factor of the deuteron 

using effective field theory. The momentum dependence of the anapole form 
factor is different from that of the matrix element of the strange currents 
in the deuteron, which may provide a way for disentangling these two com- 
peting effects when analyzing parity violating electron-deuteron scattering 
experiments. We give closed form expressions for both the form factor and 
integrals often encountered in the NN effective field theory. 
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I. INTRODUCTION 

Several open questions in nuclear parity violation Pj,0] prompt further investigation of 
parity violating observables in the simplest nuclear system. Recently, a precise measurement 
of the parity violating angular asymmetry in itp — > d'j was proposed |§ in order to tightly 
constrain the weak one-pion-nucleon coupling constant, h^ N . In addition, MIT-Bates has 
a dedicated program to measure parity violation in electron scattering on the deuteron [|J, 
which depends upon both the strangeness content of the deuteron, and on weak nonleptonic 
interactions which generate an anapole interaction. Both measurements will constrain h^ NN , 
free from finite density effects that contaminate extractions from measurements in light and 
heavy nuclei. 



*On leave from the Department of Physics, Duke University, Durham NC 27708. 
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For low momentum scattering, an effective field theory treatment fl5[-p7l of the deuteron is 
both systematic and predictive [10 .25] . The language and notation used in this paper is based 
upon refs. || and ||26|| . Here we generalize the calculation of the deuteron anapole moment 
by retaining the dependence upon the momentum transfer from the virtual photon, 



thereby determining the full anapole form factor at leading order (LO) in the effective field 
theory expansion. 

If sufficient data can be taken at Q 2 low enough for the effective field theory to be valid, 
the momentum dependence of the parity violation from electron scattering on the deuteron 
should follow, up to higher order corrections, the form we obtain. In that case, it may be 
possible to extract not only the strange axial matrix element of the nucleon, but also the 
(presently controversial) value of the parity violating pion-nucleon coupling. 



II. EFFECTIVE LAGRANGIAN 



We give here the effective Lagrangian needed for this calculation Pp6|. A discussion 



of the power counting used to determine the order at which each diagram contributes is 
presented there. Both strong and weak interactions are required for the leading order anapole 
form factor. 

The strong interactions of the pions and nucleons are described by 
Cstrong = ^TrD^D^ + [iD + N + ^iVV ■ (£D£t - £}T>£)N 

+ ^T N ^ K ° + K ^) a • BN - C^ Sl \N T P ifi N)\N T P i>0 N) + ■■■ , (2.1) 
where the pion fields are contained in a special unitary matrix, 

with / = 132 MeV. iV is the isospin doublet of nucleon fields, Mn is the nucleon mass, 
qa ~ 1-25 is the axial coupling constant, V 1 is the covariant derivative, and the ellipses 
represent operators involving more insertions of the light quark mass matrix, meson fields, 
and spatial derivatives. The isoscalar and isovector magnetic moments are k,q = 0.44 and 
K\ = 2.35 in units of nuclear magnetons. P^q is the spin-isospin projector for the 3 S\ channel, 
defined by 

1 1 

Pi,0 = —7=<72(7i Tl , ^PifiPjfi = > ( 2 - 3 ) 

V o ^ 

The coefficient Cq = —5.51 fm 2 has been determined from fits to the S-wave phase shifts 
in the 3 Si channels , renormalized at /i = m n in the power divergence subtraction scheme 
(PDS) [|. 

The leading order weak interaction Lagrange density is 

^ = i h^ NN vr + pt n + h.c. + ■ • • (2.4) 
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FIG. 1. Leading order diagrams contributing to the deuteron anapole form factor. These are 
the single nucleon contributions, written as An(\\s\) in the text. The crossed circles represent 
operators that create or annihilate two nucleons with the quantum numbers of the deuteron. The 
solid square is the weak operator with coefficient h^ NN , and the solid circle is the minimal coupling 
to the electromagnetic field. Wavy lines are photons, solid lines are nucleons, and dashed lines are 
mesons. 



where the ellipses represent higher order terms in the chiral and momentum expansion. 

The value of the parity violating pion-nucleon coupling, h^ NN , is controversial. Theo- 
retical treatments provide qualitative ranges, but experimental results are in disagreement 
|T|J^j2^j2^] . However, if we take the generally accepted size of h£l N ~ 1(T 7 |gg-i|, then 
h^ NN is responsible for the leading order parity violating effects in the deuteron. Strong 
interaction effects in the deuteron are dominated by Cq , with higher order corrections 
coming from pion exchange and contact terms multiplying operators with two derivatives 
or quark mass insertions. In the KSW power counting (outlined in II), strong interaction 
pions are treated in perturbation theory, and are sub-leading compared to the four-nucleon 
contact interactions. In the calculations of deuteron observables and interactions that have 
been computed with KSW power counting, there are no indications that the pionic con- 
tributions are larger than naively expected. Calculations at next-to-next-to leading order 
(NNLO) (for recent progress in this area see p6,3"7|) will allow more definite conclusions 
about the convergence of the theory to be made |38| . 



III. THE ANAPOLE FORM FACTOR OF THE DEUTERON 

The anapole form factor of the deuteron, Ao(|k|), describes its spin-dependent coupling 
to the electromagnetic field for arbitrary momentum transfer, defined through 



£ = iA D (\k\) -L e abc D^D%F» c , (3.1) 

N 



where k is the three momentum of the photon, F^ c is the electromagnetic field strength 
tensor, and ie a t, c D a ^ D b is the deuteron spin operator. The equations of motion for the 
electromagnetic field allow this interaction to be written in terms of local four-fermion in- 
teractions, and therefore this does not correspond to a coupling between a nucleon and an 
on-shell photon. 

The anapole form factor, Ax>(|k|), is separated into three parts: the single nucleon 
contribution An, the magnetic contribution Am, and the electric contribution Ae, with 

A D (|k|) = ^(|k|) + A M (|k|) + ^(|k|) . (3.2) 

The expressions we obtain for the form factors use methods outlined in ref. While the 



regularization and renormalization prescription used in |27| does not directly lead to gauge 
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FIG. 2. Leading order diagrams contributing to the deuteron anapole form factor. These are 
the magnetic coupling contributions, written as Aj^(\k.\) in the text. The crossed circles represent 
operators that create or annihilate two nucleons with the quantum numbers of the deuteron. The 
solid square is the weak operator with coefficient h^ N , and the open circle is the magnetic coupling 
to the electromagnetic field. Wavy lines are photons, solid lines are nucleons, and dashed lines are 
mesons. 




FIG. 3. Leading order diagrams contributing to the deuteron anapole form factor. These are 
electric coupling contributions, written as A^(|k|) in the text. The crossed circles represent oper- 
ators that create or annihilate two nucleons with the quantum numbers of the deuteron. The solid 
square is the weak operator with coefficient h^ N , and the solid circle is the minimal coupling to 
the electromagnetic field. Wavy lines are photons, solid lines are nucleons, and dashed lines are 
mesons. 



invariant results, by matching orders in a |k| 2 expansion to known dimensionally regularized 
results, we can recover the full dimensionally regulated expression. In the appendix we show 
individual integrals which are relevant for this calculation (as well as other two loop pion 
exchange calculations). 

The diagrams in fig. (JTJ) that constitute the contribution from the single nucleon factorize 
into the product of the single nucleon anapole form factor (generated by h^ N , see |39|) with 
the LO deuteron charge form factor to yield 



An(M) 



g A eh^ NN M 2 N 2 7 
TTf |k| 4 



-) tan^ 1 



mJkl 



2m„ 



tan 




(3.3) 



where 7 = y/M^B is the binding momentum of the deuteron with B the deuteron binding 
energy. 

The magnetic contribution, arising from the interaction of the electromagnetic field with 
the isovector magnetic moment of the nucleon, comes from the diagrams shown in fig. @j. 
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where Lii{z) is the dilogarithm function, given in the appendix. 

The electric contribution, arising from the minimal coupling of the electromagnetic field 
to pions and non-relativistic nucleons, comes from the diagrams shown in fig. (0). 
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(3.5) 



where p = \/|k| 2 + 4m 2 . 



IV. DISCUSSION OF RESULTS 



For illustrative purposes, we explicitly factor out the normalization of each of the form 
factors contributing to the anapole moment of the deuteron. The normalized single nucleon, 
magnetic, and electric form factors, A^, Am, and Ae respectively, are defined through 
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FIG. 4. The three contributions to the anapole form factor as a function of momentum transfer 
|k|, each normalized to unity at |k| = 0. The solid curve corresponds to the contribution from the 
single nucleon anapole form factor, yljv(|k|) (see fig. ^)). The dashed curve shows the contribution 
from the photon coupling to the isovector nucleon magnetic moment, Am{\^\) (see fig. ^)). The 
dotted curve corresponds to the electric form factor, Ae(\^\) (see fig. ^)). 
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where A M (0) = ^^(0) = A E (0) = 1. The deuteron anapole form factor, A£>(|k|), is 
normalized so that ^4^(0) = 1. (Note that while Ad = An + Am + Ae, the normalized Ad 
does not equal the sum of the individual tilded form factors.) 

Fig. (|]) shows that at LO in the effective field theory expansion the form factor asso- 
ciated with the single nucleon and magnetic contributions fall somewhat faster than the 
contribution from the electric form factor, versus increasing momentum transfer. The power 
counting for this effective field theory suggests that corrections to the LO result are 
naively at the 30% level. While it is possible that the differences between the form factors 
computed at LO are modified by higher order contributions, there is no reason for the form 
factors to be identical. 

The anapole form factor of the deuteron contributes to parity violation in electron scat- 
tering off the deuteron. A simultaneous source of parity violation in such an experiment 
is provided by any intrinsic strangeness existing in the deuteron. In a previous calculation 
|26| , we gave an expression for the leading order parity violating matrix element in electron- 
deuteron scattering in terms of the deuteron anapole moment and the strange axial matrix 
element of the nucleon. The present calculation extends this work, and increases the possi- 
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FIG. 5. A comparison between the normalized anapole form factor, Ad(\]s.\) defined in eq. (4-l)> 
and the charge form factor, expected to describe the matrix element of the strange axial current at 
LO. The solid curve is the anapole form factor while the dotted curve corresponds to the charge 
form factor Fq(\Ic\). 



bility of separating the two effects, by generalizing the moment to the momentum dependent 
form factor. To do this we need to estimate how the strange axial matrix element depends 
on momentum transfer. 

The scale for the variation of the anapole form factors is set by the deuteron binding 
momentum, 7, and the pion mass. In contrast, the variation of the form factors describing 
the matrix element of the strange axial current is set by the deuteron binding momentum, 
and the kaon mass or some higher mass scale. Therefore, we expect the contribution to 
parity violating electron-deuteron scattering from the matrix element of the axial strange 
quark current to become more important at higher momentum transfers. The fall of the 
strange quark form factor with momentum should be driven by the small deuteron binding 
energy; that is, by the nucleon recoil effects. This gives the LO strange deuteron form factor 
the same structure as the LO deuteron charge form factor, 

Fc(|k|) = ^tan-(M) . (4 .2) 

The parity violating amplitude in electron-deuteron scattering is then, at LO in the effective 
field theory, 

A<ri{\k\) = 10- 7 (s.4 g AS F c (\k\) - 0.21 A D (\k\) j ^ e M e 7i e e* a e h . (4.3) 

The two contributions to the anapole form factor are shown in fig. (^). A more precise 
estimate of the form factors, and hence a more precise estimate of h^ N through a compar- 
ison with electron-deuteron scattering data, would require a higher order calculation of the 
anapole form factor in addition to more insight into possible strange quark effects (while 
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current measurements of g& s are consistent with zero HO- 42], so is one of the measurements 
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An estimate of the size of the higher order corrections entering at NLO and 



higher can be made by examining the size of the corrections to the charge form factor, as 
presented in [10|- Compact, analytic expressions for the LO form factors have been obtained 
in this work, and it is likely that similarly simple expressions can be obtained at NLO. 

It is appropriate to comment on the recent work by Khriplovich and Korkin who 
determine the anapole moment of the deuteron in the zero-range potential limit. They find 
(using their revised expression for the magnetic moment term) a contribution from the elec- 
tric and magnetic interactions (not including the single nucleon contribution) proportional 
to 



.4 



(KK) 
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"V +7 

(m* + 2 7 ) 2 



(4.4) 



The expression we obtained in |26f does not reproduce their result, except in the limiting 
case where 7 — > 0. It is clear that their expression, shown in eq. ( [4.4| ), does not have the 
correct behavior in the chiral limit for fixed 7 (the deuteron anapole moment would diverge, 
driven by the single nucleon contribution), as discussed in |[26|| . The deuteron anapole 
moment cannot have any l/m n dependence in the chiral limit when 7 is held fixed, due 
to the vanishing of the strong one-pion-deuteron coupling. Our result, shown in eq. ( |4.1| ). 
correctly reproduces this necessary behavior. 



V. CONCLUSION 

In this work we have presented an analytic expression for the anapole form factor at 
LO in the effective field theory expansion. The dominant contribution to this form factor 
is from pion physics, through the weak pion-nucleon coupling h^ NN . The pion mass and 
deuteron binding momentum are the scales that drive its variation. In contrast, pion effects 
will not dominate the matrix element of the strange axial current. We find that while the 
deuteron binding momentum provides most of the momentum variation of both form factors, 
the anapole form factor falls fasters than the matrix element of the strange axial current. 
This has important implications for the determination of both g& s and h^ NN from parity 
violating electron-deuteron scattering. Encouraged by the success of effective field theory 
for this parity violating process, and the relatively simple closed form expressions we obtain, 
we are optimistic about higher order calculations of this form factor and the possibilities of 
extracting h^ NN from electron scattering experiments. 

This work is supported in part by the U.S. Dept. of Energy under Grants No. DE-FG03- 
97ER4014 and DE-FG02-96ER40945, and NSF grant number 9870475. 
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VI. APPENDIX OF INTEGRALS 



The analytic expressions we found for the anapole form factors come from integrals which 
are relevant for other calculations. In this appendix we give expressions for the general form 
of some of the integrals encountered in evaluating the Feynman diagrams for the anapole 
calculation. We use a combination of dimensional regularization and the position space 
representation (as advocated by Binger |27j) to find expressions for the integrals. The 
divergent contribution and well defined finite piece will be given in terms of 

r = r[4-n] +iog(7r) + i , (6.1) 

where n is the number of space-time dimensions. To be consistent with the PDS procedure, 
fi/2 and not \x is used for the dimensional regularization mass-scale in the definition of the 
integrals. In the expressions that follow, the magnitude of the three-momentum transfer 
is given by k = |k|. All integrals have been power-series expanded about n = 4, so terms 
of order (n — 4) 1 or higher are not shown. Where needed to distinguish one integral from 
another within a class (same denominators in the integrand), a subscript is used to indicate 
the naive degree of ultraviolet divergence of the diagram. 



A. 1-Loop, 2-Propagators 

We find that two-loop integrals often reduce down to some simple functions, dilogarithms, 
and the following integral (which also appears in one loop diagrams): 

r d 3 q 1 1 1 _ x I k \ 

A ^ k ) = ] ( 27 r)3[( q -Ik) 2 + a 2 ][q 2 + & 2 ] = 2^fc tan [W+Qj ' C ' 



B. 2-Loops, 3- Propagators 

A frequently encountered integral in the theory with perturbative pions |8||| is 
B (a,b,c; fj,;k) evaluated at k = (an integral that is ultra-violet divergent as n — > 4) 
where 



B (a,b,c;n; k) 
2 



8-2n , d n~l ^n-lj 



(27r) n - 1 (27T)"" 1 [(q - |k) 2 + a 2 } [(q + l) 2 + b 2 } [l 2 + c 2 ] 



(6.3) 



It can be computed at k = using dimensional regularization: 
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(6.4) 
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where M = a + b + c. The methods of p7 are used to find the finite fc-dependent portion, 
which matches onto the previous result to give the full expression for arbitrary k: 



-Bo (a, b, c; /i; k) 



32tt 2 



r-iot 



M 2 + 
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k 
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2M 



+ 2 



(6.5) 



This reduces to eq. ( |6.4j) in the limit /c — > 0, as required. 
A related integral is 
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2, 



q ■ k 
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(6.6) 



Making the replacement q- k — > — (q— |k) 2 + q 2 + |k 2 leaves a more divergent integral (the 
q 2 term) than found in the B integral, but a judicious combination of terms allows explicit 
evaluation of this integral using the methods of |[27|| . We find 
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(6.7) 



C. 2-Loops, 4-Propagators 



In graphs where an external current couples to a nucleon line and a perturbative pion 
is exchanged between nucleons, we encounter integrals with four propagators and two loop 
momentum integrations. Such integrals are expressed in terms of simple functions and 
dilogarithmic functions. The function P(X, Y\ k) is defined to be 
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where the dilogarithmic function Li 2 (z) is conventionally defined 
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The simplest integral involving four propagators at two loops does not have any loop 
momentum dependence in the numerator, and is finite at n — 4. It is expressed in terms of 
functions defined previously, 



7_ 2 (a, b, c, d; k) 
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(6.10) 



An integral with a more complicated momentum dependence, but without Lorentz struc- 
ture in the numerator of the integrand, can be written in terms of i_ 2 , 
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We will use the following integrals in subsequent expressions. 

d 3 q dH q-k 1 1 
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The first integral with a Lorentz index, but finite at n — 4, is 
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(6.14) 



where we have used the known external momentum dependence to rewrite the 1 integrand 
in terms of the q variable. After completing squares in the numerator we obtain 
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The following integrals have two indices in the numerator of the integrand. These again 
are naively divergent and must be treated in n dimensional space-time. The divergences all 
enter through integrals already defined, B , B ± , J_ 2 , and I . 
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8-2n , rf n-l q ^n-q q .j 



1 

~~ 2 
and 



c — d 



r a 11 q gT x l q 

J (27T)"- 1 (27T)"- 1 [(q- |k) 2 + a 2 ] [q 2 + 6 2 ] [(q + l) 2 + c 2 ] [l 2 + d 2 ] 

^ A(a, b; k) — B (a, c, d; /j,; k) + (d 2 + 5 2 — c 2 ) 7-2 (o, c, d; fc) 



(6.17) 



J\~o(a, c, d; //; fc) 

1 //A 8 " 2 ™ /■ d n - x q d"- 1 ! (q-k)(l-k) 1 



k 2 V2 
1 

2P 



(27T)™- 1 (27T)™- 1 [(q - ±k) 2 + a 2 ] [q 2 + b 2 } [(q + l) 2 + c 2 ] [l 2 + d 2 ] 
d — ^ a_b_ + ^ 2 _ fl2 _ b- kyj - Si (a, c, d; /x; fc) 



47T 



A: 2 



— (d + 6 — c ) I B (b : c, d; /x; 0) — Bo(a, c, d; k) + (b —a — — )i_ 2 (a, fe, c, d; k) 



d — c 
An 



a 2 + ^ 

A(0, b; 0) - A(a, 6; fc) ^ (A(a, 0; fc) - A(a, 6; k)) 

k 2 



+B (0, c, d; n; 0) - -B (o, c, d; k) — (a + — ) J_ 2 (a, 0, c, d; fc) 



_(rf 2 + 6 2 _ c ^ 



— (-B o (0, c, d; 0) - B Q (b, c, d; 0)) - J_ 2 (a, 6, c, d; fc) 



a 2 _|_ 

1 (7_ 2 (a, 0, c, d; fc) - i_ 2 (a, 6, c, d; fc)) 



b 2 



(6.18) 
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It is necessary to keep the n-dependence of the coefficients in the tensor structure of eq.( |6.16| ) 
explicit, since factors of n — 4 combine with divergences in the K Q and K to produce finite 
contributions. These are necessary to recover a gauge invariant result for the weak one- 
photon matrix element in the deuteron, and hence the anapole moment. 
Finally we present another two index integral 



L 



S-2n 



d n - L q d n ~ l \ 

(27T)™- 1 (27T)"- 1 

q* q J 



1 



n 



[(q - ±k) 2 + a 2 } [(q + ±k) 2 + b 2 } [(q + l) 2 + c 2 ] [l 2 + d 2 } 
L - Lq) S<> + — ^ ((n - 1)L - L ) 1 U 



n 



k 2 



(6.19) 



where 

L (a, b, c, d; \x\ k) 



8-2n 



q^ 



(27T)"" 1 (27T)"" 1 [(q - ±k) 2 + a 2 ] [(q + ±k) 2 + 6 2 ] [(q + l) 2 + c 2 ] [l 2 + d 2 } 



k 2 

B (b, c, d; fi; k) + B (a, c, d; //; fc) — {b 2 + a 2 + — ) J_ 2 (a, 6, c, d; fc) 



(6.20) 



and 



L (a, 6, c, d; k) 

1 (fi\ 8 - 2n f d n ~ l 
~ k 2 V2 



q d n -H 



(q • k) 2 



(27T)™- 1 (27T)"- 1 [(q - |k) 2 + a 2 ] [(q + ±k) 2 + b 2 } [(q + l) 2 + c 2 ] [l 2 + d 2 } 



1 

2T 2 



l k 2 



a* 



\9 



! k 2 



U 2 



Io(a, y — + y + — , c, d; fi; k) + !„(&, d — + — + — , c, d; fi; k) 



(6.21) 



We have shown the set of integrals needed to calculate the anapole form factor of the 
deuteron. We see that the diagrams involved, even with finite momentum transfer from 
external currents, are computible in closed form. These results will be useful for many other 
two-loop computations. 
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